Power dissipation for systems with junctions of multiple quantum wires 
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We study power dissipation for systems of multiple quantum wires meeting at a junction, in 
terms of a current splitting matrix (M) describing the junction. We present a unified framework 
for studying dissipation for wires with either interacting electrons (i.e., Tomonaga-Luttinger liquid 
wires with Fermi liquid leads) or non-interacting electrons. We show that for a given matrix M, the 
eigenvalues of M T M characterize the dissipation, and the eigenvectors identify the combinations of 
bias voltages which need to be applied to the different wires in order to maximize the dissipation 
associated with the junction. We use our analysis to propose and study some microscopic models of 
a dissipative junction which employ the edge states of a quantum Hall liquid. These models realize 
some specific forms of the M-matrix whose entries depends on the tunneling amplitudes between 
the different edges. 

PACS numbers: 73.23.-b, 73.40.Gk, 71.10.Pm 



I. INTRODUCTION 

One-dimensional systems of strongly correlated elec- 
trons have been studied extensively for several years both 
experimentally, in the form of quantum wires and carbon 
nanotubes, and theoretically pHy]. Junctions of several 
quantum wires have also been studied in recent years 
since they can now be experimentally realized in carbon 
nanotubes The existing studies of junctions of 

quantum wires, which are usually modeled as Tomonaga- 
Luttinger liquids (TLL), have mainly looked at their low- 
temperature fix ed p oints and the corresponding conduc- 
tance matrices |r3 - l24| . Many of these studies have fo- 
cused on situations in which there is no power dissipation 
in the system. The aim of our work will be to include 
dissipation in the discussion. For simplicity, we will con- 
sider only spinless electrons and will restrict ourselves to 
the zero frequency limit (DC) in our work. 

A motivation for studying dissipation is as follows. In 
Ref. [l6l a two-parameter description of a junction of 
three TLLs has been discussed. The dissipationless fixed 
points were shown to lie on the circumference of a circle, 
while the interior of the circle corresponds to dissipative 
junction. In this context, the center of the circle which 
corresponds to the current splitting matrix with all its 
elements equal to 1/3 is of particular interest as it cor- 
responds to the maximum possible dissipation allowed 
by constraint of current conservation. It therefore seems 
useful to understand dissipative junctions in a general 
way and to study whether any of the points inside the 
circle correspond to fixed points of some renormalization 
group (RG) equations. 

The plan of this paper is as follows. In Sec. II, we intro- 
duce the idea of a current splitting matrix M at a junction 
of N wires and discuss the cases of both non-interacting 
and interacting electrons. In Sec. Ill, we obtain an ex- 
pression for the power dissipated in terms of this ma- 



trix. The measure of the degree of dissipation is then 
defined in terms of the eigenvalues of M T M. For two- 
wire and three-wire junctions, we write down the most 
general form of the M-matrix allowed by current con- 
servation, thus providing a complete parametrization of 
the dissipation at the junction. In general, an M-matrix 
which respects current conservation can have both posi- 
tive and negative matrix elements. As we will show, an 
S-matrix describing non-interacting electrons scattering 
at the junction can be related to a matrix M all of whose 
elements are positive; this relation will follow from the 
assumption that there are no phase correlations between 
electrons coming from different reservoirs which lie far 
away from the junction. But when M has negative el- 
ements, such a relation does not exist and the matrix 
then necessarily corresponds to a system of interacting 
electrons. 

In Sec. IV, we introduce a simple model involving three 
patches of the edge states of a quantum Hall liquid with 
filling fraction v. The patches are taken to be mutu- 
ally coupled to each other by local electron tunnelings 
between three distinct points lying on the three patches 
with amplitudes cr^ , where i,j denote the patch index. 
Then a parametrization of the M-matrix is obtained in 
terms of the conductance amplitudes ct^-. In this way 
we obtain interesting dissipationless matrices in the lim- 
its Uij — > and oo respectively. These matrices were 
shown to represent dual fixed points in the theory of a 
junction of TLL wires in Ref. using more involved 
calculations. In Sec. V, we introduce a more complex 
model of a junction of three quantum wires in which the 
junction consists of a ring-shaped region with edges of 
its own. Once again, the matrix M of the entire system 
can be found in terms of the coupling of each wire to the 
ring and the tunneling amplitudes across the two edges of 
the ring. Even though such a geometry is complicated, 
from an experimental point of view it allows for easier 
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tunability as far as realizing various types of M matri- 
ces is concerned. In Sec. VI, we make some concluding 
remarks. 

II. THE CURRENT SPLITTING MATRIX 

A junction is a meeting point of N wires each of which 
has an incoming and an outgoing mode. Physically, if the 
junction is made of a material like a carbon nanotube, 
then the incoming and outgoing modes (which carry cur- 
rents) belonging to a single wire are not separated in 
space. But if these are quantum wires made out of the 
edge states of a quantum Hall liquid [25[ , then the incom- 
ing and outgoing modes are spatially separated. In the 
following discussion, we will consider a junction of several 
quantum wires, each with two spatially separated chiral 
current carrying edges, one incoming and one outgoing. 
Each chiral mode (incoming or outgoing) is labeled by 
an index i which runs from 1 to N and is parametrized 
by a coordinate x. We will take x to run from to oo; 
the point x = will be common to all the wires and 
will denote the junction. The outgoing currents in the 
system are related to the incoming currents by a current 
splitting matrix M given by 

Joi = £ M y Jij- (1) 

3 

Current conservation at the junction therefore implies 
that each column of M must add up to 1. Let us also 
assume that the incoming current on wire i is propor- 
tional to the applied bias voltage Vn, with the constant 
of proportionality being the same for all wires. Then if 
all the wires have the same bias voltage, the net outgo- 
ing current Joi — Jn on each wire i must vanish which 
implies that each row of M must also add up to 1 . 

For non-interacting electrons, the junction can be de- 
scribed in terms of a scattering matrix § which provides a 
linear relation between the incoming and outgoing elec- 
tron fields at the junction. Namely, the incoming and 
outgoing electron fields, ipn(x,t) and if)oi{x,t), are re- 
lated at all times t as 

iMM) = § v iMM)- ( 2 ) 

3 

Current conservation implies that § must be an N x N 
unitary matrix. Any deviation from the linear boundary 
condition for the electron fields at the junction will imply 
the existence of local inter-electron interactions at the 
junction even if the electrons in the bulk of the wires are 
left non-interacting. The scattering matrix description 
can also be used to describe electrons which are weakly 
interacting in the bulk of the wire, by trea ting the effects 
due to interactions perturbatively [lj, HE G3] • 

Given a scattering matrix S for non-interacting elec- 
trons, we will now see how the elements of the current 
splitting matrix M can be found. The incoming and out- 
going currents Jn and Joi in wire i are proportional 



to \ipn\ 2 and | "0oi 1 2 respectively. Eq. ([2]) implies that 

Hoi} 2 = J2jk^ij^ ik ^h^' Ik - We now assume tnat tnere 

are no phase correlations between the incoming electrons 
on different wires j and k since they are coming from 
different reservoirs whose distances from the junction 
are taken to be much larger than the phase coherence 
length; the absence of such phase correlations is crucial 
for the validity of the Landauer-Biittikcr theory of elec- 
tronic transport in mesoscopic systems [26j j . Hence terms 
like ipljtpik can be set equal to zero if j ^ k. We thus 
obtain IV'Oil 2 = J2j l^iiPlV'/jl 2 - This is 01 the same form 
as in Eq. ([T]) if we identify My = ISy) 2 - 

On the other hand, if we have strongly interacting 
electrons in one dimension, then it is natural to use 
bosonization. The electrons in the wire are then ex- 
pressed in terms of free bosonic excitations described 
by TLL theory, and the fixed point theory of the junc- 
tion can be described in terms of a current splitting ma- 
trix M which is obtained by imposing a linear boundary 
condition on the incoming and outgoing bosonic fields 
at the junction Q3, QJ, [ll [llQ]. One can use free 
bosonic fields to describe either non-interacting or in- 
teracting electrons in the bulk of the one-dimensional 
wires depending on whether the Luttinger parameter g 
is equal to or not equal to 1. Note that even when 
we have g = 1 (non-interacting electrons) in the bulk 
of the wire, within the bosonization approach the cur- 
rent splitting matrix M representing a linear relation 
between the incoming and outgoing boson fields at the 
junction corresponds to the presence of non-zero inter- 
electron interaction at the junction. This is because the 
boson fields are related to the corresponding electron 
fields by a non-linear bosonization identity, ipi/o( x ) — 

(l/V^Tra) F 1/O e l ' t, " o{x) [H-@|, where Vj/oO) are the in- 
coming and outgoing electron fields, <pj iq are the in- 
coming and outgoing chiral bosonic fields, and Fjio are 
the corresponding Klein factors. Thus there is a subtle 
difference between using a scattering matrix S for non- 
interacting electrons and a current splitting matrix M 
in bosonized TLL theory for electrons which are non- 
interacting (g = 1) in the bulk of the wire. In the lat- 
ter case the M-matrix description of the junction corre- 
sponds to an interacting theory of electrons where the 
interaction is localized at the junction. Now, if the in- 
coming and outgoing boson fields are linearly related 
to each other at the junction, i.e., if 4>cn(x — 0,i) = 
J2j (x = 0,t), then M must be a real and orthogo- 

nal field splitting matrix in order that both the incoming 
and outgoing bosonic fields satisfy the canonical com- 
mutation relations [l8|, H3, HI] . Such a description of the 
junction given by an orthogonal M represents fixed points 
of the junction as was shown in Ref. [TH The current at 
any point of wire i is given by —(l/2n)d(/)/dt. Hence 
we note that the above condition at the junction implies 
that the outgoing and incoming currents also satisfy Eq. 
([T]), i.e., the field splitting M matrix can be taken to be 
the same as the current splitting matrix. Hence in the 
bosonic formalism, M must be an N x N real and orthog- 
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onal matrix each of whose rows and columns add up to 

i m. 

As we will see later, the orthogonality condition on 
the M matrix also renders it dissipationless irrespective 
of its origin, i.e., this is true for both a junction of non- 
interacting electrons described by an § matrix or a junc- 
tion of TLL wires described by a field splitting matrix M. 
Hence, to include dissipation in the analysis we have to 
relax the condition of orthogonality on M. This also im- 
plies that the formalism of bosonization cannot be used 
directly since a non-orthogonal M does not allow the 
bosonic commutation relations to be satisfied. However, 
M must continue to be real since it relates incoming and 
outgoing currents which are all real, and each of its rows 
and columns must add up to 1 as argued before. Hence 
the main emphasis of this section is on the fact that var- 
ious situations comprising of either non-interacting elec- 
trons or interacting electrons, where the interaction is 
either localized at the junction or extended all over the 
wire, can be described just in terms of a current splitting- 
matrix (Eq. (|T]|). We will see in the following section that 
this information is enough to characterize the dissipation 
associated with the junction. 



III. DISSIPATION 

We will now consider the specific case involving edge 
states in a quantum Hall liquid for discussing dissipation 
in a junction. In such systems, currents flow only along 
the edges as all states in the bulk are localized; such 
edge modes are chiral in nature and can be described 
by theories of chiral bosons [29[. In the linear response 
regime, if a voltage V is applied to an Ohmic contact 
which is assumed to be perfectly coupled to the edge, 
then for filling fraction v, the current J injected into 
the edge from that contact is given by J = GV , where 
G = ve 2 /h is the conductance. 

Now let us derive an expression for the power dissi- 
pated by such a system which is governed by a current 
splitting matrix M. The power Pd dissipated near the 
junction is given by the difference of the total incoming 
and outgoing powers [3(| HJ , 

i 

= ^ JJ (I M T M) Jj, (3) 

where we have introduced a matrix notation in the sec- 
ond line of Eq. ((3|) , with Jj being a column made up of 
the incoming currents Ju, and I being the N x N identity 
matrix. On physical grounds, the power dissipated near 
the junction cannot be negative. This implies that the 
eigenvalues Xi of M T M must necessarily lie in the range 
[0,1]. If the incoming current Jj is proportional to an 
eigenvector of M T M with eigenvalue Xi, the power dissi- 
pated will be proportional to 1 — A*. The set of values 



of 1 — A, therefore provides a measure of the amount of 
dissipation associated with a system characterized by the 
current splitting matrix M. 

We emphasize here that Eq. ([3]) describes the power 
dissipated in the region close to the junction, and not in 
the leads which are assumed to be far away from the junc- 
tion. For instance, for a two-wire junction with perfect 
transmission of the currents, i.e., for a matrix M given 
by Mn = M 2 2 = and Mi 2 = M21 = 1, the expression 
in Eq. ([3]) vanishes; however, we know that dissipation 
occurs in the leads because the outgoing electrons eventu- 
ally equilibrate to the chemical potential there, leading 
to a contact resistance of e 2 /h (for spinless electrons). 
Thus, M T M = / only means that there is no dissipa- 
tion associated with the junction, although dissipation 
can still occur in the leads. 

Since each row and column of M adds up to unity, 
both M and M T must have one eigenvalue equal to 1, 
the corresponding eigenvector being given by a column 
all of whose entries are equal to each other. This col- 
umn is therefore an eigenvector of M T M with eigenvalue 
equal to 1 which corresponds to a situation where the 
bias voltages Vn (or incoming currents Ju) on all the 
wires are equal, and no power is dissipated. Also note 
that the power dissipated vanishes for all possible values 
of the incoming currents if M is orthogonal. On the other 
hand, the dissipated power is maximized if all the eigen- 
values of M T M are equal to except for one eigenvalue 
which is necessarily equal to 1. This occurs when all the 
entries of M are equal to 1/N. Hence if we think of a 
situation where the dissipation happens at the junction 
and not in the leads, the entire incoming power will be 
converted to heat at the junction and the outgoing power 
will vanish. 

In general, a current splitting matrix M corresponding 
to an iV-wire junction has (N — l) 2 independent param- 
eters. This is because the first (N ~ 1) x (A — 1) block of 
M can have arbitrary entries while the entries of the last 
row and column of M are then fixed by the conditions 
that each row and column must add up to 1. For N = 2, 
we need only one parameter and the matrix is given by 

«=(iVr). w 

where a must lie in the range [0, 1] to ensure that the 
dissipated power is always non-negative. No power will 
be dissipated if a = or 1 (i.e., M is orthogonal), while 
the maximum power can be dissipated if a = 1/2. It is 
interesting to note that this one parameter family of M- 
matrices can be obtained from the following electronic S- 
matrix describing scattering of non- interacting electrons, 

s = / VT ^). (5) 

The case of maximum dissipation, i.e, a = 1/2 for the 
non-interacting electrons case also corresponds to ex- 
tremal shot noise [32j as is expected. 
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For N = 3, we require four parameters to specify M in 
general as we can see below 

(a b 1 — a — b \ 

c d 1-c-d . (6) 

1 — a — c 1 — b — d a + b — c — d— 1 / 

(The ranges of the parameters a — d are fixed by the con- 
dition that the power dissipated must be non-negative; 
hence we will not specify these ranges here). If we de- 
mand that no power be dissipated, i.e., that M be or- 
thogonal, then we only need to specify one parameter as 
will be discussed below. Note that the three-wire case is 
quite different from the two- wire case discussed earlier. 
In the three- wire case, it was possible for M to have some 
negative elements without violating current conservation 
and non- negativity of the dissipated power, in sharp con- 
trast to the two- wire case. To get a better feel for this, let 
us consider a one-parameter family of M-matrices which 
corresponds to a highly symmetric junction given by 

/ a (l-a)/2 (l-a)/2 \ 
M = (1 -a)/2 a (l-a)/2 • ( 7 ) 
\ (1 -a)/2 (l-a)/2 a / 

This matrix corresponds to a situation in which the re- 
flected current in each wire and the transmitted currents 
from one wire to the other two are the same for all the 
wires. Using the condition of non- negativity of the net 
dissipated power, we can show that the parameter a must 
lie between —1/3 and 1. For a = —1/3 and 1, M is or- 
thogonal and is therefore dissipationless. For all values 
of a lying between —1/3 and the diagonal elements of 
M are negative. It is easy to see that M-matrices with 
negative entries cannot be obtained from any unitary Si- 
matrix, i.e., cannot be obtained from any non-interacting 
electron theory. Hence such M-matrices necessarily cor- 
respond to situations in which the inter-electron interac- 
tion strength is non-zero. We emphasize that such cur- 
rent splitting matrices only exist for a junction of three 
or more wires and are absent for the two- wire case. 

Next, let us consider a situation where the power asso- 
ciated with the incoming current is set to unity in units of 
1/(2G). Then the three-element column given by Jj can 
be identified with a unit vector in three dimensions which 
can be parametrized as (sin 9 cos <p, sin 9 sin 0, cos 9). The 
maximum dissipation occurs for the M-matrix which 
has all its elements equal to 1/N. This corresponds to 
a = 1/3 in Eq. Q. Using Eq. © for the case of a = 1/3, 
we find that the power associated with the outgoing cur- 
rents is given by 

P — — [cosQ + sin6(cos(j) + sm^)] 2 . (8) 
o 

Note that P a is bounded by [0, 1] and is symmetric under 
9 — > 7r — 9 and <f> —></> + n. To visualize the expression 
in Eq. flSJ), we present a contour plot of P a as a function 
of <t> and 9 in Fig. [TJ An interesting point to note in 
the figure is the existence of a line of points on which 




FIG. 1: (Color online) Contour plot of P given in Eq. (JHJ) as 
a function of <j> and 6. 



the outgoing power is zero and hence the power dissi- 
pation is maximum. This implies that there is a fam- 
ily of bias voltage or incoming current configurations for 
which the power dissipated is maximum. On the other 
hand, there are two points at which the outgoing power 
is unity which correspond to zero power dissipation. To 
understand these patterns, we recall that the direction 
of maximum power dissipation in the space of incoming 
current vectors corresponds to the two distinct eigenvec- 
tors of the M T M with zero eigenvalue. For any value of 
a in Eq. 0, an orthonormal set of eigenvectors of M T M 
is given by Vi = (1, 1, 1)/V3, V 2 = (1, -1, l)/y/2, V 3 = 
(1,1, — 2)/v6 and the corresponding eigenvalues are 1, 
(1 — 3a) 2 /4 and (1 — 3a) 2 /4. Note that these eigenvectors 
are independent of the parameter a. This is so because of 
the symmetric form of M matrix. Hence, for this entire 
family of M-matrices (Eq. (0), the eigenvectors (i.e., the 
combinations of incoming currents) which give the max- 
imum power dissipation are independent of a. Second, 
the eigenvalue which is different from unity is a quadratic 
function of a which is zero for a = 1/3 (maximum dis- 
sipation) and unity for a = —1/3 and a — 1 (both cor- 
responding to zero dissipation). The line of maximum 
dissipation appearing in Fig. [T] corresponds to an in- 
coming current column J/ which is a linear combination 
of the two degenerate eigenvectors corresponding to the 
zero eigenvalue given by J/ = cosS V 2 + sind V 3 , where 
S lies in the interval [0, 27r]. The existence of such a line 
of maximum dissipation is encouraging from an experi- 
mental point of view since this implies that we only need 
to vary a single parameter in an experiment to encounter 
the point of maximum dissipation. The two points in Fig. 
[1] which have P a = 1 (zero dissipation) correspond to the 
eigenvector Vi. There are two such points because we 
get zero dissipation if the incoming current is prepared 
either in the direction of this eigenvector or opposite to it; 
in Fig. [Q these points lie at (0,0) = (tt/4, cos _1 (1/\/3)) 
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and (57r/4, tt — cos' - X (l/V3)). To conclude, we see that a 
study of the eigenvectors of M T M can lead to a complete 
understanding of dissipation in a junction as a function 
of the bias voltages applied in the various wires. 



IV. A THREE- WIRE MODEL WITH 
DISSIPATION 

In this section, we develop a microscopic model for 
a three-wire system with a dissipative junction. A 
schematic picture of the system is presented in Fig. [5] 
The currents and voltages on each wire will be assumed 
to be governed by J = GV (where G — ue 2 /h) on all the 
incoming and outgoing chiral wires. (The symbols Vi in 
the figure denote the incoming voltages which drive the 
incoming currents; the outgoing currents and voltages are 
then determined by the Vi and the matrix M which will 
be derived below). The junction region consists of three 
points a, 6, c, one point lying on each of the three wires as 
shown in Fig. [2] Electrons can tunnel between any two 
of these points, say, i and j. If the tunneling amplitude is 
denoted by £y , the corresponding tunneling conductance 
crijG will be proportional to |£y| ■ Here we have intro- 
duced the quantity G so that <xy is dimensionless. The 
conductances satisfy = (Tji > 0. If the voltages at the 
two points are given by Vi and Vj, the current flowing 
from i to j will be given by (JijG(Vi — Vj). We will now 
see that this model gives rise to a current splitting matrix 
M which is generally dissipative. 




FIG. 2: (Color online) Picture of a three-wire model with 
tunneling conductances atj between points i,j which can take 
values a, b, c. V% denote the incoming voltages. 

In our analysis, we will work directly with the currents 
without introducing any fermionic or bosonic fields. To 
derive the matrix M, we have to determine the outgoing 
currents (Joi, J02, Jos) m terms of the incoming cur- 
rents (J/i, J12, J/3)- The incoming and outgoing currents 
(and therefore voltages) will generally change discontin- 
uously at the three junction points. The corresponding 



incoming and outgoing voltages are obtained by divid- 
ing the currents by G. We assume that the voltages at 
each of the three points of the junction are given by the 
mean values of the corresponding incoming and outgoing 
voltages. Namely, 



V 



1 



( V n + V 0i ) 



1 

2G 



( Jn + J 0l ) (9) 



for i = 1,2,3. 

The mean value assumption made in Eq. (JSJ can be 
justified as follows. We can begin with a model in which 
the tunneling region is not a point but has a finite length 
I, and there is a tunneling conductance per unit length 
given by dij [27l l33j . Tunneling will then occur from 
every point lying in the tunneling region in wire i to the 
corresponding point lying in wire j. We then find that 
the current Ji on wire i changes smoothly from Jn to 
Joi as we go from one end of the tunneling region to 
the other. Hence the voltage Vi(xi) = Ji(xi)/G will also 
change smoothly, where i runs over 1,2,3, and Xi runs 
over the tunneling region from to The current Ji(xi) 
can be obtained by solving equations of continuity given 

by[27|,H 



r) J 

^ = - a tJ [JiM-Jjixj)], (10) 



We can solve these equations to obtain the dependence 
of the current Ji(xi) and voltage Vi(xi) — Ji(xi)/G on 
the coordinates x%. If we now take the limit I — > with 
Idij = Uij being held fixed, we recover the earlier model 
of tunneling between three points, with the voltages at 
the three points being given by Eq. ©. 

We now return to our original model and write down 
equations of continuity for the currents at the three tun- 
neling points, 



Joi 



■hi — 



G 



(Vi - Vj), (11) 



for i = 1, 2, 3. Using Eq. ©, we can solve for the Joi in 
terms of the Jn. This enables us to obtain the matrix M 
which relates the two sets of currents. We find that 



ill 



1 + <J23 - (l/4)5 2 



Si 



012 



+ (3/4)S 2 ' 
(l/2)5 2 



1 + Si + (3/4)5 2 



(12) 



where Si = 012 + 023 + 031 and S 2 = 012023 + 023031 + 
031 012- All the other entries of M can be found by sym- 
metry. Note that each row and column of M adds up to 
1 as desired. In addition, M being a symmetric matrix is 
a special feature of this specific model. 

We can show in general that the expression for power 
dissipation given in Eq. Q agrees with the sum of 
the powers dissipated by the three tunneling processes. 
Namely, if we substitute the expression for M given in Eq. 
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(fT2| in Eq. ([3]) , and compare that with the expression for 
the power dissipated (= I X V) by the three tunnelings, 
namely, 

G[a 12 (V 1 -V 2 ) 2 +a 23 {V 2 -V 3 ) 2 +a 31 (V 3 -V 1 )% (13) 

(where Vi — Vj appears in Eq. (jlljl ). we find that the two 
agree for all values of the incoming currents Jn. 

In the special case that o\ 2 = a 23 — a 3 \ = a, the 
expression for M simplifies to 



1 



1 + 3cr/2 




(14) 



which is of the form given in Eq. (|7|)- Three particular 
values of this M are worth noting, namely, 



for (7 = 0, 




1/3 1/3 1/3 \ 

1/3 1/3 1/3 for a = 2/3, 
1/3 1/3 1/3 J 

-1/3 2/3 2/3 
2/3 -1/3 2/3 
2/3 2/3 -1/3 



for a = oo. (15) 



I- matrix given in Eq. (|14p. one of the eigen- 
T M is equal to 1, while the other two are 
3cr)] 2 . There is no power dis- 



For the 
values of 

equal to [(2 - 3<x)/(2 
sipation (M is orthogonal) if a = or oo, while there 
is maximum power dissipation if a = 2/3. On physical 
grounds it is natural to expect that there is no dissipa- 
tion if er = 0. But the dissipation also turns out to to be 
zero for a = oo which is somewhat surprising. This can 
be traced back to the analysis done in Ref. [l6| for the 
junction of three TLL wires. In that paper, the authors 
started with a situation where there is a perfectly reflect- 
ing (disconnected) junction of three TLL wires effectively 
described by a M-matrix corresponding to a = in Eq. 
(1141) . and then switched on electron tunneling operators 
between each pair of wires such that the amplitudes of all 
the three tunneling operators are equal. Using the tech- 
nique of bosonization, they then established that as the 
strengths of all the tunneling operators go to infinity un- 
der an RG flow, the system is described by the M-matrix 
which is obtained by taking the a = oo limit in Eq. (1141) . 
In their analysis, the M-matrices corresponding to both 
(7 = and (7 = oo are fixed points of the theory, and they 
are connected to each other by a duality transformation. 
These statements make our model seem quite attractive 
because even though it is rather simple, it manages to 
capture the essential non-trivial physics related to dissi- 
pation in a three-wire junction without getting into the 
technicalities of bosonization. 

If the various edges shown in Fig. [5] are the edges of a 
quantum Hall system, the tunneling operators will satisfy 



some RG equations. Depending on the filling fraction v 
and the location of the quantum Hall liquid with respect 
to the edges [llj], the tunneling operators will be either 
irrelevant or relevant, and the corresponding tunneling 
conductances will then flow to or oo respectively. This 
implies that the RG fixed point for M will be given by 
either the first matrix or the last matrix in Eq. (|15l) . 
Thus, the second matrix in (1151) which corresponds to 
maximum dissipation does not appear to be a fixed point 
of the model. 



V. A MORE COMPLEX THREE- WIRE MODEL 

We now consider another model for a dissipative junc- 
tion of three wires. This model consists of a ring shaped 
region (with two chiral edges) and three external wires 
(each with two chiral edges: incoming and outgoing) 
which connect to the ring at three different points. All 
the edges carry currents and can be modeled by TLLs. 
Further, each of the external wires can have different bias 
voltages which determine the incoming currents imping- 
ing on the ring region. Along the ring, the co-propagating 
currents can tunnel between the two edges. For simplicity 
we will assume equal tunneling amplitudes at all points. 




FIG. 3: (Color online) Picture of a three- wire model where the 
external wires connect to a ring which has two co-propagating 
edges with interedge tunneling. Vi denote the incoming volt- 
ages. 

Fig. [3] gives a schematic picture of the model we have 
in mind. Each of the external wires is made of a TLL (or 
a single edge of a fractional quantum Hall system) and 
has an outgoing and an incoming chiral edge Ioi and In, 
where i labels the wire. The ring is also made of a TLL 
(or a single mode fractional quantum Hall edge) and has 
two co-propagating modes, one on the outer edge (Joi) 
and the other on the inner edge (Jn). 
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At each of the 'point' junctions where an external wire 
meets the ring, we have three incoming modes and three 
outgoing modes (marked by arrows in Fig. [3]). Such a 
three- wire junction can be described by a orthogonal 3x3 
current splitting matrix (M) whose rows and columns 
add up to 1. The orthogonality implies that (i) the junc- 
tion relates the outgoing bosonic fields to the incoming 
bosonic fields in a way which preserves the chiral commu- 
tation relations of the fields, and (ii) the junction is dis- 
sipationless. For a three-wire charge-conserving and dis- 
sipationless junction, the matrix M can be parametrized 
by a single continuous parameter 9 [llj 0, [28| , and it can 
be classified into two classes for which (a) detMi = 1, 
and (b) detMb = — 1. These two classes are expressed as 



(16) 



In Eq. ([Til), a = (1 + 2cos6>)/3, b = (1 - cos6> + 
v / 3sin6')/3, and c = (1 — cos9 — V3sin0)/3. In the 
Mi class, 9 = corresponds to the disconnected N fixed 
point, 9 = 7r to the Dp fixed point, and 9 = ±27r/3 to 
the chiral fixed points x± m the notation of Ref. UM 

Each of the three 'point' junctions in Fig. [3] is char- 
acterized by a dissipationless current splitting matrix M. 
For simplicity we will now assume all the three junctions 
have the same M with the same orientation. We will also 
assume all the 'point' junction matrices to be identical 
and of the Mi type. Next, we will allow tunneling be- 
tween the inner and outer edges of the ring, which can be 
thought of 'classically' as a resistor connecting the inner 
and outer current carrying wires. A more microscopic 
model of such a dissipative tunneling is given in Refs. [27| 
and K33l The main result is that at the ends of each tun- 
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neling region (of length L), the currents on the outgoing 
edges are a linear combination of the incoming currents 
and can be written as 



Joi(L) 
Jo2(L) 



1-t t 
t 1-t 



Jn(0) 
J/ 2 (0) 



(17) 



The parameter t can be expressed in terms of the micro- 
scopic tunneling conductance as 



t 



(1 _ e -2Lah/(ve 2 )^ 



(18) 



for the case of co-propagating edges, where a is the tun- 
neling conductance per unit length between the two edges 
of the ring (27| . 

We note again that the 'point' junction matrices con- 
necting external wires to the ring are dissipationless. The 
only source of dissipation in our model is therefore the 
interedge tunneling between the co-propagating modes 
propagating on the ring. Now, starting from a given dis- 
sipationless current splitting matrix Mi at each 'point' 
junction and a given interedge tunneling parameter t, we 
can solve for the three outgoing and six interedge cur- 
rents in terms of the three incoming currents. We then 
find that the M-matrix of the system which relates the 
outgoing currents to the incoming currents is of the cyclic 
form 



(19) 





where d, e and / are given by 



30i 2 - 48< + 27 + (60t 2 - 84t + 42) cos 9 + (l8t 2 - 30t + 12) cos(2(9) 
42t 2 - 48t + 33 + (28t 2 - 68t + 34) cos 9 + (38t 2 - A6t + 14) cos(2(9) ' 

12t 2 + 6t- (24< 2 - 12t + 6) cos 9 + (l2t 2 - 18t + 6) cos(26>) 
42< 2 - 48t + 33 + (28i 2 - 68t + 34) cos 9 + (38i 2 - 46i + 14) cos(26») ' 
-6t + 6 - (8i 2 - At + 2) cos 9 + (8t 2 + 2t - 4) cos(26>) 
f ~ 42< 2 - 48t + 33 + (28i 2 - 68t + 34) cos 9 + (38i 2 - 46£ + 14) cos(26») ' 

If we take all the matrices at the three 'point' junctions to be identical and of the M2 type, we again find that the 
M-matrix of the complete system is of the cyclic form given in Eq. p^|) . although the expressions for d,e,f are 
different from those given in Eq. (f2"U|). For the M-matrix given in Eqs. (|19H20p . one of the eigenvalues of M T M is 
equal to 1 (non-dissipative), while the other two (degenerate and dissipative) are given by 

78t 2 - 84t + 33 + (28t 2 - 68t + 34) cos 9 + (2t 2 - lOt + 14) cos(26>) 
A ~ 42i 2 - 48t + 33 + (28i 2 - 68t + 34) cos 9 + (38t 2 - 46i + 14) cos(26>) ' ^ 



r 



A contour plot of A in the t — 9 plane is presented in Fig. |4j Since A is symmetric under 9 — > —9, we have only 
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plotted 9 from to 7r in the figure. We see that A = 1 
(no dissipation) if either 9 = 0, it or t = 0. It was shown 
in Ref. [27| that an RG flow takes the variable La to 
either or oo, depending on the value of the interaction 
parameter of the TLLs which constitute the two edges of 
the ring. Hence the fixed-point values of the parameter 
t are and 1/2 according to Eq. (fT5|> . For t — > 0, the 
eigenvalue A goes to 1 for any value of 9 and we therefore 
get a dissipationless M-matrix. But for t — > 1/2, we find 
that 

_ 21 + 14cosl9 + 19cos(26») 

39 + 14cos0 + cos(20) ^ ' 

This is equal to 1 for 9 — 0, tt, and is not equal to for any 
value of 9. According to Fig. @J the point of maximum 
dissipation (A = 0) lies at (t,9/w) ~ (0.419,0.583), and 
not at t = or 1/2, and it is therefore not a fixed point 
of this model. 



' ^^1 
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FIG. 4: (Color online) Contour plot of A given in Eq. (|2T|) as 
a function of t and 6. 



VI. DISCUSSION AND SUMMARY 

In this work we introduced a scheme to quantify dis- 
sipation for a A^-wire junction for both non-interacting 
electrons and TLL wires. The quantification is achieved 
in terms of a real current splitting matrix M. The dissi- 
pated power can be parametrized by the non-zero eigen- 



values of I — M T M and hence there is no dissipation if M 
is orthogonal since / — M T M is then equal to 0. We have 
shown that if an eigenvalue of I — M T M is equal to 1, 
the corresponding eigenvector determines a combination 
of the applied bias voltages for which the input power is 
completely dissipated at the junction. For a three-wire 
junction, the matrix M with all entries equal to 1/3 has 
a doubly degenerate eigenvalue equal to 1. Hence any 
linear combination of the two eigenvectors corresponds 
to a combination of bias voltages which will lead to com- 
plete dissipation at the junction. This implies that the 
bias voltage combination which corresponds to maximum 
dissipation is not a unique point in the allowed parame- 
ter space but forms a one-parameter family of points as 
discussed in Sec. III. This fact makes it more likely to be 
accessible in an experimental situation. 

We presented two microscopic models of dissipation 
for a three- wire system, one involving tunneling between 
three points (Sec. IV) and the other involving tunneling 
between three pairs of edges lying on a ring (Sec. V). 
The model in Sec. IV leads to a symmetric M-matrix 
depending on three parameters a^ , while the model in 
Sec. V leads, for a particular choice of current splitting 
matrices at the 'point' junctions, to a cyclic M-matrix 
depending on two parameters t, 9. 

For both models, we have briefly discussed the RG 
flows of the various parameters. For the model in Sec. 
IV, the RG flow takes the system to one of two fixed 
points, both of which correspond to dissipationless M- 
matrices. For the model in Sec. V, the RG flow again 
takes the system to one of two fixed points, one of which 
gives a dissipationless matrix while the other is generally 
dissipative (except for the special cases 6 = 0,%). In all 
cases, we find that the matrix corresponding to maximum 
dissipation (i.e., all elements of M being equal to 1/3) is 
not a fixed point of the RG equations. Hence within the 
models we have studied, it appears that there is nothing 
special about the maximally dissipative M-matrix from 
an RG point of view. 
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